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Abstract 

The diffusive relaxation of a colloidal fluid adsorbed in a porous medium depends on many factors, 
including the concentration and composition of the adsorbed colloidal fluid, the average structure 
of the porous matrix, and the nature of the colloid-colloid and colloid-substrate interactions. A 
simple manner to describe these effects is to model the porous medium as a set of spherical particles 
fixed in space at random positions with prescribed statistical structural properties. Within this 
model one may describe the relaxation of concentration fluctuations of the adsorbed fluid by simply 
setting to zero the short-time mobility of one species (the porous matrix) in a theory of the dynamics 
of equilibrium colloidal mixtures, or by extending such dynamic theory to explicitly consider the 
porous matrix as a random external field, as recently done in the framework of mode coupling 
theory [V. Krakoviack, Phys. Rev. Lett. 94, 065703 (2005)]. Here we consider the first approach 
and employ the self-consistent generalized Langevin equation (SCGLE) theory of the dynamics of 
equilibrium colloidal mixtures, to describe the dynamics of the mobile component. We focus on 
the short- and intermediate-time regimes, which we compare with Brownian dynamics simulations 
involving a binary mixture with screened Coulomb interactions for two models of the average static 
structure of the matrix: a porous matrix constructed by quenching configurations of an equilibrium 
mixture in which both species were first equilibrated together, and a pre-existing matrix with 
prescribed average structure, in which we later add the mobile species. Wc conclude that in both 
cases, if the correct static structure factors are provided as input, the SCGLE theory correctly 
predicts the main features of the dynamics of the permeating fiuid. 

PACS numbers: 05.40.-a, 82.70.Dd, 02.70.-c 
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Many relevant systems and processes in industry and nature involve the diffusion of 
colloidal dispersions through porous media The most relevant issues involve either 
the equilibrium and phase behavior j^, or the transport and dynamic properties 3], of 
the permeating colloidal liquid. Thus, one would like to understand, from a fundamental 
perspective, how these properties depend on factors such as the porosity and morphology 
of the matrix, the nature of the interaction forces (among the particles and with the porous 
matrix), or the concentration of the colloidal dispersion. Some of these issues require the 
development or extension of the fundamental description of the dynamic behavior of bulk 
colloidal systems 0, 0] to the case in which these systems permeate a porous matrix. 

The fundamental study of the phenomena above rely on simplified models. Thus, a porous 
medium is sometimes modeled as a simple geometry (planar slit, cylindrical pore, etc.) to 
describe local phenomena, whereas random arrays of locally regular pores incorporate the 
intrinsic randomness of most natural or synthetic porous materials 1|. One may adopt, 
instead, a simplified model of a random porous medium, namely, a matrix of spherical 
particles with random but fixed positions. This matrix is permeated by a colloidal liquid, 
whose dynamics we wish to understand. Such model systems have been employed to describe 
mostly equilibrium structural properties ^], although simple model experimental realizations 
of this system have been prepared [?], in which one could also measure the dynamic properties 
of the mobile species. The interpretation of such measurements requires sound theoretical 
schemes to describe the dynamics of the permeating dispersion. One possible approach uses 
available theories of the dynamic properties of bulk colloidal mixtures j^, Iq, 1^ in which the 
mobility of one of the species is artificially set equal to zero. Another possibility is to first 
reformulate these theories to explicitly consider the porous matrix as a random external field 



In this letter we demonstrate that the first of these approaches suffices to correctly 



predict the main features of the dynamics of the permeating fluid, provided the correct 
average static structure of the matrix and of the adsorbed fluid is available. 

This conclusion is based on the use of the multi-component self-consistent generalized 
Langevin equation (SCGLE) theory of colloid dynamics [9], to describe the relaxation of 
concentration fluctuations of the mobile component in the model porous matrix. By setting 
the free-diffusion coefficient of one species to zero, this theory is readily adapted to the de- 
scription of the dynamics of the model system above, thus allowing the numerical calculation 
of the partial intermediate scattering function F{k,t) of the mobile species. As a concrete 
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and illustrative application, here we report the theoretical predictions for a binary colloidal 
mixture of particles interacting through screened Coulomb potentials in which one of the two 
species plays the role of the porous matrix. The theoretical results for this specific system 
are compared with the corresponding results of a Brownian dynamics (BD) simulation on a 
model system consisting of = N1 + N2 Brownian particles in a volume V, with A''i(= niV) 
particles of species 1 and A^2(= "^2^) particles of species 2 interacting through direct forces, 
but not through hydrodynamic interactions. 

We carry out two kinds of computer experiments that differ in the manner we generate 
the structure of the porous matrix. In the first kind, we let the particles of both species 
';o undergo Brownian motion, according to the conventional Brownian dynamics algorithm 



12] with the same free- diffusion coefficient = D2 = until equilibrium is reached. At 
this point, we artificially arrest the motion of the particles of species 2 by setting D2 = 
at an arbitrary configuration. In the second kind of experiments, a pre-existing matrix 
is formed in the absence of the mobile species, by choosing the arrested configurations 
according to a prescribed distribution, afterward "pouring" the mobile particles into this 
matrix of obstacles. The prescribed average structure of the matrix that we consider below 
corresponds to the structure of an equilibrium mono- component fluid of species 2. In both 
cases, after choosing a particular configuration of the matrix, we let the mobile species 
equilibrate in the external field of the fixed particles at that particular frozen configuration, 
and then proceed to the calculation of the dynamic properties of interest. Since these 
properties depend on the specific configuration of the fixed particles, and in order to average 
out this dependence, the results presented here correspond to an average over more than 
50 different configurations of the porous matrix. In both cases we also record the radial 
distribution functions between the two species, to be employed as the static input required 
by the SCGLE theory. In our illustrative application, the direct interactions are described 
by a hard-sphere plus a strong repulsive screened Coulomb interaction. For simplicity, we 
assume that both species have the same hard-sphere diameter a, so that the potential (in 
units of the thermal energy fc^T = (3~^) is given by Puais^r) = +00 for r < a, and for r > cr 
by 



puai3{r) = JK^Kp . (1) 

V ^jg- 
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One may think of the parameter \fK^ as proportional to the electric charge of species 
a, and the parameter z = kg y& the dimensionless inverse Debye length k jsl- 

The relevant dynamic information of an equilibrium zz-component colloidal suspension is 
contained in the v^v matrix F{k,t) whose elements are the partial intermediate scattering 
functions Fap{k,t) = (^^(k, t)r;,^(— k', 0)) where ^^(k, t) = I]^°iexp[ik- ri{t)]/ y/Na, with 
rj(t) being the position of particle i of species a at time t. The initial value Fa[3{k,0) 
is the partial static structure factor Sap{k) j^, Q. In our simulation experiment, we are 
interested in the dynamic properties of the mobile species, represented by F{k, t) = Fii{k, t). 
This dynamic property will be theoretically calculated applying the multi-com^nent self- 
consistent generalized Langevin equation (SCGLE) theory of colloid dynamics [9,] with the 
particular condition u = 2 and = 0. 

n n 

The SCGLE theory, explained in more detail in Refs. [9|, 110[, is summarized by a self- 
consistent system of equations for the uxu matrices F{k, t) and F^^\k, t) (the latter defined 
as F'^^{k^t) = 6a[3 ^exp [ik • AR(")(t)]^, where AR'^"^(i(:) is the displacement of any of the 
Na particles of species a over a time t, and 6ai3 is Kronecker's delta function. Written in 
matrix form and in Laplace space, and omitting the explicit /c-dependence, the self-consistent 
system of equations reads 

Fiz) = {z + {I+ [AC{z)]\)-'eDS-'y' S, (2) 

and 

= {z + {I+[AC{z)]\)-'k'Dy\ (3) 

where 5* is the matrix of partial static structure factors, D and X{k) are diagonal matrices 
given by Da/s = Sa/sD^ and Xa/sik) = Sai3[l + {k / k[°'^Y]^^ , where fc^"^ is the location of the 
first minimum (following the main peak) of Saa{k)- AC*(t) is a diagonal matrix with its 
diagonal element A(*(t) given by 

AC(t) = d'kk''[F^'\t)U[cV^Fit)S-'V^hU, (4) 

where the elements of the /c-dependent matrices h and c are the Fourier transforms haj3{k) 
and Cajjik) of the Ornstein-Zernike total and direct correlation functions, respectively. Thus, 
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h and c are related to 5* by S* = / + ^Jnh^fn = [/ — \fnc^Jn\ ^, with the matrix ^Jn defined 

as = 5a(3^/n^. 

The control parameters of our system are the interaction parameters z, Ki and K2, 
and the volume fractions (pi and 02 (with (p^ = irriaCr^/Q). Here we fix the value of the 
screening parameter to z = 0.15, and start by considering an equi-molar mixture with 
01 = 02 = 2.2 X 10^*^. In our first simulation experiment we start with the simplest 
case, namely, a mono-disperse suspension of = A^i + A''2 identical particles interacting 
with the same pair potential {Ki = K2 = 100), which execute Brownian motion. After 
thermalization, we stop the motion of half of them, and let the other half constitute the 
mobile species. Two additional experiments of the same kind were performed for systems 
with the same parameters as above, but varying the coupling parameters Ki and K2. Thus, 
the second experiment corresponds to a more interacting system, Ki = K2 = 500, and 
the third to an asymmetric mixture such that the matrix is formed by the more strongly 
interacting particles, K2 = 500 and Ki = 100. These three experiments belong to the first 
kind referred to above, i.e., they involve a matrix whose average static structure is identical 
to the partial static structure factor S22{k) of an equilibrium mixture of both species. In this 
kind of experiments, 6*22 (^) and the other static structural properties may be determined 
during the initial equilibration stage, before arresting the motion of the matrix. As an 
example, in Fig. [Ha we present the various radial distribution functions gajsi^) simulated in 
this manner in the third of these experiments. 

We also performed parallel experiments of the second kind, involving a pre-existing matrix 
with prescribed average structure. The prescribed structure we chose corresponds to the 
equilibrium structure of a mono- component system containing only species 2. Thus, we first 
let the A^2 particles equilibrate, and then freeze an arbitrary configuration in which we then 
place the A'^i particles of the other species. The static and the dynamic properties involving 
the mobile species are then simulated after the fiuid of species 1 equilibrates in the external 
field of the matrix in that particular configuration. The results are then averaged over a 
sufficient number (~50) of configurations of the matrix. In Fig. [Hb we present the resulting 
fi'o/3(^) corresponding to the third experiment of this second kind, involving a system with 
the same parameters as in Fig. [Ha. 

These simulated structural properties are now employed as the static input needed by the 
SCGLE theory, thus avoiding the use of liquid state approximations [l3|. The predictions of 
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FIG. 1: Brownian dynamics simulated radial distribution functions gapir) of a colloidal fluid 
(species 1) diffusing through a porous matrix formed by a second species of fixed particles, in- 
teracting with the repulsive Yukawa potential with fixed screening parameter z = 0.15, volume 
fractions 4>i = 4>2 = 2.2 x 10~^ and repulsion strength parameters Ki = 100 and K2 = 500 with 
the porous matrix generated in the presence (a), and in the absence (b), of the mobile species. 

the SCGLE theory for the dynamic properties of the Brownian fluid immersed in the porous 
matrix are presented in Fig. [2] for the three experiments of the first kind and for only the 
last experiment of the second kind, namely, that involving the asymmetric mixture with 
K2 = 500 and Ki = 100. These results are compared with the corresponding BD results 
for the normalized time-dependent diffusion coefficient D*{t) =< (Ar(t))^ > /{QD^t) and 
for the collective intermediate scattering function F{k,t) = Fii{k,t) of the mobile species 
at fixed k and for the times t = 0, to, and lOto, with to = a'^ / . Notice that for t = lOto, 
D*(t) has relaxed from its initial value of 1 towards, and close to, its asymptotic value, 
characteristic of the long-time regime. Thus, the illustrative data in this figure covers the 
so-called short- and intermediate-time regimes, most easily accessible by Brownian dynamics 
simulations or video-microscopy experiments. 

From information such as that summarized in Fig. [2] we may highlight the following. 
First, the description of the SCGLE theory for symmetric systems (illustrated by the first 
two rows) is highly accurate in the time-regimes illustrated in the figure. Second, the 
comparisons in Fig. [2] have essentially the same quality as the corresponding comparisons 
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FIG. 2: Time-dependent self diffusion coefficient D*{t) (left column) and partial intermediate 
scattering function F{k, t) for t = 0, to, and lOto (right column) of the diffusive species permeating 
the porous matrix, interacting with the repulsive Yukawa potential with fixed screening parameter 
z = 0.15 and volume fractions </>! = i;^>2 = 2.2 x 10~^, but with parameters Ki and K2 given by 
Ki = K2 = 100 (first row), Ki = K2 = 500 (second row) and Ki = 100 and K2 = 500 (third and 
fourth rows). The symbols represent Brownian dynamics results and the solid lines are the SCGLE 
theoretical predictions. 

involving fully thermalized mixtures in which both species diffuse jol. Third, the theoretical 
predictions for F{k, t) in asymmetric systems for both kinds of experiments (last two rows) 
have similar levels of quantitative accuracy; the largest discrepancies with the simulation 
data are illustrated by the results of the fourth experiment, without being particularly severe. 
Clearly, any improvement will involve either reformulating the intrinsic approximations of 
the SCGLE theory or developing an extension similar to that carried out by Krakoviack ll| 
in the framework of mode coupling theory (MCT) At best, however, such improvements 
will only add to the quantitative accuracy of the present results and will be of the magnitude 
illustrated in Fig. [21 

Thus, the comparisons above indicate that the SCGLE theory, devised to describe the 
dynamics of equilibrium colloidal mixtures, provides a useful approach to the dynamics 
of a mono-disperse suspension permeating a porous medium formed by a random array 
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of other colloidal particles. This approach may now be applied to explore other interesting 
phenomena such as, for example, transitions of dynamic arrest of colloidal mixtures in porous 
media Q|. In fact, the present theory, complemented by adequate liquid state approximations 
for the static structure jl3^, may be used as a fully theoretical first-principles approach to 
qualitatively scan other regions of the state space to locate interesting dynamic phenomena 
for which simulations or experiments would be difficult to carry out, or are not yet available. 
In doing this we only loose quantitative precision, but no qualitative accuracy. In fact, one of 
the objectives of adapting the SCGLE theory to the description of the dynamics of colloidal 
dispersions adsorbed in porous media is the description of dynamic arrest phenomena in 
these systems. The analysis of its accuracy in the short- and intermediate times presented 
here is an important step in the process of assessing its overall reliability. We mention, 
however, that the SCGLE theory has been successfully employed to describe dynamic arrest 
in bulk mono-disperse Q, IgI and multi-component l3] colloidal systems. In fact, we 
can communicate that the SCGLE theory presented here also predicts dynamic arrest in 
the systems considered above upon increasing the coupling parameters Ki and/or K2 or 
by varying other control parameters of the system. In this manner, we can outline the full 
dynamic arrest phase diagram of this model system. In particular we can report that for 
systems involving only har d sp here interactions, and under the specific conditions studied 
by Krakoviack using MCT , we obtain essentially the same dynamic arrest scenario that 
he derives from MCT. Our approach, however, opens the possibility to study more complex 
situations in a simpler manner and from an alternative perspective to that provided by 
MCT. The details of these specific advances, however, will be reported separately. 
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